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ABSTRACT 


In  this  report  we  investigate  the  effect  of  real¬ 
istic  propagation  paths  on  the  visual  and  spectral 
amplitudes  of  fundamental  mode  Rayleigh  waves.  We  show 
that  the  Harkrider  amplitude  response  factor  of  a 
layered  media  should  relate  directly  to  a  station 
correction;  it  varies  by  a  factor  of  nearly  three  be¬ 
tween  oceanic  and  continental  sites.  Attenuation  due 
to  effective  Q  losses  may  cause  scatter,  of  over  one- 
half  magnitude  unit,  in  teleseismic  amplitudes  but  is 
unimportant  at  regional  distances.  Differences  in 
recorded  amplitude  caused  by  dispersion  over  oceanic 
and  continental  paths  are  sufficient  to  warrant  the 
application  of  path  corrections  to  visually  measured 
amplitudes.  These  differences  can  be  avoided  by  comput¬ 
ing  spectra  of  the  signals,  or  suppressed  somewhat  by 
using  the  results  of  stationary-phase  approximations. 
More  uniform  propagation  path  effects  at  periods 
longer  than  20  seconds  suggests  that  use  of  40  or  50 
second  waves  would  provide  better  seismic  discrimina¬ 
tion  capability  and  better  yield  estimation. 


INTRODUCTION 

In  the  estimation  of  surface-wave  magnitude,  the 
raw  data  are  Rayleigh  wave  signals  at  several  stations. 
These  time  series  must  be  condensed  into  single  values; 
in  addition  to  choosing  a  particular  cycle  of  motion  on 
the  time  record  or  a  particular  frequency  point  in  the 
spectrum  to  give  this  value,  certain  corrections  for 
the  propagation  paths  are  necessary.  The  choice  of 
where  to  measure  is  not  totally  objective  and  systematic 
in  seismological  practice,  and  the  corrections  for  the 
propagation  paths  are  not  fully  defined.  Observed 
Rayleigh  waves  at  different  azimuths  but  the  same 
distance  from  an  event  can  vary  considerably  in  visual 
character  and  by  large  factors  of  amplitude  at  any 
particular  frequency. 

A  recorded  surface-wave  signal  can  be  represented 
schematically  as  the  end  result  of  the  process; 


where  each  step  involves  a  convolution  in  the  time 
domain.  Because  of  the  more  easily  describable  and 
predictable  nature  of  long-period  surface  waves,  as 
opposed  to  short-period  body  waves,  we  are  able  to 
investigate  the  effects  of  the  source  and  the  path  on 
Mg  determinations  in  a  relatively  complete  manner.  The 
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source  effects  are  represented  by  radiation  patterns  of 
energy  from  equivalent  force  systems  and  have  been 
discussed  by  von  Seggern  (1969),  who  showed  how  Mg  of  an 
earthquake  is  dependent  upon  depth  of  focus,  azimuth 
from  the  source  to  the  recording  station(s),  and  the 
orientation  of  the  double-couple  forces.  It  was  con¬ 
cluded  that  the  distribution  of  stations  and  the  source 
geometry  could  cause  some  earthquakes  to  appear  like 
explosions  when  the  Mg  vs  discriminant  was  applied. 

In  this  report  we  investigate  the  effects  of 
travel  path  on  estimates  of  Mg  by  examining  amplitude 
response  factors,  attentuation  losses,  and  dispersion 
characteristics  of  different  earth  structures.  Specif¬ 
ically,  we  desire  to  know  whether  a  factor  dependent  on 
propagation  path  can  be  removed  from  the  observed  ampli¬ 
tudes  of  Rayleigh  waves  so  that  the  variance  of  Ms 
determinations  at  several  stations  for  a  single  event 
can  be  reduced  and  so  that  a  single  Mg  determination 
would  be  more  closely  related  to  yield  or  seismic  moment. 
Since  it  has  already  been  shown  that  the  variance  of  m^ 
is  almost  twice  that  of  Mg  (von  Seggern,  1970;  Ericsson, 
1971)  ,  we  believe  it  is  better  to  use  surface  waves 
(or  long-period  body  waves)  to  estimate  those  parameters. 
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SEPARATION  OF  SURFACE  WAVE  SIGNALS 
INTO  SOURCE,  PATH  AND  RECEIVER  FACTORS 


Our  representation  follows  that  of  Sat®  (i960)  and 
is  a  straightforward  convolution  of  all  the  time-domain 
elements  which  affect  the  surface  wave  signal,  or 
alternatively  a  multiplication  of  frequency-domain 
elements.  For  Rayleigh  waves  propagating  in  a  layered 
earth,  the  frequency-domain  expression  of  a  recorded 
signal  is 

W(u>)  *  S(u)  L(u)  B(u)  I  (u) 


where 


S  (u)  -  source  spectrum 
L  (u>)  *  layer  response 
B  (u>)  ■  attenuation  spectrum 
I  (u>)  *  instrument  response 

We  consider  each  of  these  factors  separately  in  the 
following  sections.  We  will  also  examine  dispersion, 
which  is  actually  part  of  the  layer  response  but  which 
deserves  to  be  considered  separately  because  it 
strongly  affects  the  shape  of  recorded  waveforms. 

Source  spectrum 

We  consider  as  a  source  a  simple  point  force, 
whose  time  variation  has  Fourier  transform  F  (u) , 


(1) 
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directed  downward  on  the  surface  of  the  earth  along  the 
positive  axis;  the  source  spectrum  is  simply 

S(w)  -  F(w) 

To  model  more  realistic  sources  we  would  use 

S(w)  ■  P(w)a^  (2a) 

for  explosions  where  P  (u)  is  the  Fourier  transform  of 
the  pressure  time  history  at  the  equivalent  cavity 
radius  £  and  we  might  use 

S(u)  -  F(w)d  (2b) 

for  earthquakes  where  F  (u>)  is  the  Fourier  transform  of 
the  force  time  history  of  the  simple  components  of  the 
couple  or  double  couple  with  moment  arm(s)  of  length  d. 

Note  that  both  these  realistic  source  functions  have 
dimensions  of  force-length-time  as  opposed  simply  to 
force-time  for  a  simple  point  force;  therefore,  for 
the  layer  response  to  them,  a  factor  with  units  of 
length"*  must  appear,  multiplying  the  layer  response  to 
the  single  point  force.  This  factor  takes  the  form  of 
the  wavenumber  k  times  some  dimensionless  expression 
for  the  radiation  pattern. 

Layer  response 

Layer  response  functions  to  various  source  models 
for  realistic  earth  structures  have  been  described  by 
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Harkrider  (1964)  and  Ben-Menahem  and  Harkrider  (1964). 
The  vertical  component  response  to  a  single  downward 
point  force  at  or  very  near  the  surface  of  the  earth  is 
given  by  (Harkrider,  equation  (103)): 

L(u>)  =  AR(w)H^2)(kr)e“i7T/2 

where  AR  is  the  Rayleigh  amplitude  factor  of  Harkrider 
and  Il£  ^is  the  zero' th-order  Hankel  function  of  the 
second  kind.  Beyond  about  a  wavelength  from  the  source, 
the  expression  can  be  well  approximated  by 

Hg2^(kr)  *  ( 2 / tt k r )  ^2exp[-i(kr-ir/4)  ] 

so  that 

L  (id)  =  AR(a))k‘1/2r“1/2expI-i(kr+n/4)] 

where  the  constants  have  been  absorped  in  the  amplitude 
factor  Ar  (tu).  For  the  explosive  source  within  2  or  3  km 
of  the  surface,  the  layer  response  can  be  closely 
approximated  by 

L(oi)  =  -eQ((D)AR((D)k1/^2r’^//2exp[-i(kr-3TT/4)  ] 

where  eQ  is  the  ellipticity  of  the  Rayleigh  wave  at  the 
surface,  taken  to  be  negative.  Note  the  180°  phase 
difference  between  the  response  to  a  downward  point 
force  and  the  response  to  a  spherical  compressional 
source.  Note  also  that  the  additional  k  factor  preserves 


(3a) 
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the  correct  units.  For  the  earthquake  source*  the 
corresponding  medium  response  is 

L(w)  »  AR(w)kly,2r'1/,2exp[-i(kr-3Tr/4)]x(w) 

where  x  (w)  is  Ben-Menahem  and  Harkrider' s  radiation- 
pattern  factor  which  was  studied  by  von  Seggern  (1969) 
for  its  effects  on  magnitude  estimates. 

We  are  interested  in  the  structure-dependent 
factors  of  the  layer  response  to  an  explosive  source 
in  equation  (3a):  these  are  AR,  eQ,  and  k1  2.  Outputs 
for  numerous  structures  from  HARKRIDER,  a  program  which 
calculates  Rayleigh-wave  phase  and  group  velocities  as 
well  as  AR,  were  available  at  the  SDL.  For  these 
different  structures  which  include  oceanic,  shield,  ^ 
and  tectonic  types  we  found  that  no  value  of  eQ  or  k 
at  a  period  of  20  seconds  differed  from  another  by  more 
than  101.  These  two  factors  can  be  considered  to  be 
constant  for  all  structures,  and  a  magnitude  (Mg) 
correction  for  structure  need  not  involve  them.  For 
periods  longer  than  20  seconds, eQ  and  k1/2  varied  less 
among  the  structures,  and  for  shorter  periods  they 
varied  more.  The  factor  AR  at  20  seconds,  however, 
varied  between  2.82  x  10-4  for  a  deep-ocean  structure 
and  8.47  x  10-4  for  a  continental  structure  with  a 
very  thick  crust.  Harkrider  and  Anderson  (1966)  give 
the  formal  expression  for  AR  as 

ar  -H  WH 1 


(3b) 
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where  C  and  U  are  Rayleigh-wave  phase  and  group  velocity 
respectively,  p  is  density,  z  is  the  depth,  and  the 
squared  quantities  are  ratios  of  horizontal  and  vertical 
particle  velocities  at  depth  z  to  vertical  particle 
velocity  at  the  surface.  The  integral  portion,  as  well 
as  C  and  U,  vary  with  structure  and  mode.  Using  all 
the  available  HARKRIDR  outputs  at  the  SDL,  we  show 
fundamental-mode  AR  (at  20  seconds  period)  versus  CU 
(at  20  seconds  period  also)  in  Figure  1,  The  predicted 
inverse  dependence  is  evident,  with  upper  points 
representing  thick-crust  continental  structures  and 
the  lower  two  points  representing  deep-ocean  structures. 
A  simple  straight  line  (visually  adjusted)  seems  to  fit 
the  data  adequately,  since  no  point  deviates  from  it 
by  more  than  105  on  the  vertical  scale.  Plots  of  AR 
versus  U  in  Figure  2  and  AR  versus  C  in  Figure  3  show 
reasonable  straight-line  approximation  also,  with 
scatter  somewhat  greater  than  for  AR  versus  CU.  Thus, 
knowing  either  group  or  phase  velocity  of  the  medium, 
or  preferably  both,  one  could  estimate  the  factor  AR 
directly  within  about  105.  Otherwise,  one  must  invert 
the  velocity  data  to  obtain  a  structure  and  then  com¬ 
pute  Ar  for  the  structure. 

A  further  consideration  in  applying  AR  corrections 
is  that  the  propagation  path  is  rarely  a  laterally 
homogeneous  structure.  For  instance,  what  AR  is  to  be 
applied  to  a  signal  recorded  on  a  continent,  having 
originated  in  another  continent,  but  having  traversed 
an  oceanic  path  for  most  of  the  cpiccntral  distance? 

In  the  past  Rayleigh  wave  propagation  through  gradual 
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changes  in  structure  at  the  continental  boundaries  have 
been  studied  with  models  and  approximated  by  wedges  or 
even  vertical  discontinuities  in  theoretical  solutions. 
Actual  observations  such  as  that  by  McGarr  (1969)  show 
a  signal  being  amplified  by  a  factor  of  2  or  3  in 
crossing  from  oceanic  to  continental  structure  and 
deamplified  by  the  same  factor  for  the  opposite  direc¬ 
tion.  We  point  out  that  the  AR  for  oceanic  structures 
is  half  as  great  as  the  average  for  continental  struc¬ 
tures  (Figure  1),  and  this  factor  alone  may  be  a  simple 
means  of  predicting  amplitude  changes  from  one  structure 
to  another.  Moreover,  if  the  value  of  AR  at  a  point  on 
a  gradually  changing  structure  controls  the  amplitude 
at  that  point,  it  may  be  necessary  to  know  only  AR  at 
the  receiving  station  in  order  to  apply  the  AR  correc¬ 
tion  to  magnitude.  Only  the  eventual  solution  of  the 
Hay lcigh-wavc  propagation  problem  in  changing  crust- 
mantle  structures  by  finite-element  schemes  can  provide 
the  correct  answers. 

We  regard  the  factor  AR  as  primarily  a  station 
correction.  If,  however,  the  structure  in  the  source 
area  immediately  grades  into  quite  dissimilar  types, 
such  as  at  a  continental  boundary,  partition  of  the 
source  energy  into  the  various  modes  nay  vary  signifi¬ 
cantly  with  azimuth.  Considering  now  a  source  area  homo¬ 
geneous  to  a  radius  of  at  least  a  wavelength,  such 
partition  of  energy  would  he  the  same  at  all  azimuths, 
and  observed  amplitudes  of  signals  which  are  recorded 
on  different  structures  could  he  equalized  by  dividing 
then  by  this  amplitude  factor  AR  appropriate  to  the 


recording  site.  Using  the  extreme  continental  values  of 
AR,  from  approximately  4  x  10"4  to  8  x  10"4,  the  largest 
correction  fo  Mg  would  then  amount  to  the  logarithm  of 
their  ratio,  0.3,  when  no  oceanic  sites  were  involved. 
However,  in  practice,  these  extreme  structures  are 
seldom  found,  and  we  expect  typical  M  corrections 
based  on  AR  to  he  much  less  than  0.3  when  no  oceanic 
sites  arc  involved. 

Attcntuation 

Bruno  (1962)  showed  that  the  attenuation  factor 
for  surf  >cc  waves  is  given  by 

B(u)  «  exp|-ur/2Q(w)U(w) | 

where  Q  is  the  dimensionless  quality  factor  and  U  is 
the  group  velocity.  Here  Q  is  understood  to  be  the 
"effective"  Q  over  the  path  in  question.  Thus,  measured 
Q  for  Kayleigh  waves  of  20  second  period  includes 
losses  due  to  causes  other  than  actual  absorption;  e.g., 
node  conversion  and  the  various  "optical"  processes 
(reflection,  refraction,  diffraction,  scattering), 

Measured  values  for  20  second  Kayleigh  waves 
are  not  numerous;  Tsai  and  Aki  (I960)  found  that  Q 
averaged  approximately  ?oo  for  many  paths  from  the 
I’arkfield  earihquake  of  1906,  Tryggvason  (1961)  found 
to  he  over  lOiiU  for  Asian  and  European  paths, 

Marshal  and  i!atpenter  (I9t4i)  found  it  to  he  about  If 0 
for  the  JSovaya  2cm I ya  to  £orth  America  path.  In  these 
studies,  the  waves  often  traversed  major  structural 
boundaries,  von  ?*eggern  (197(1)  found  an  average  ()  of 


134  in  the  period  range  o£  XO  to  16  seconds  for  paths 
from  NTS  explosions  to  stations  around  the  Western 
United  States;  this  low  value  probably  reflects  the 
tectonic  nature  of  the  earth's  crust  and  upper  mantle 
there.  From  Figure  7  one  can  estimate  the  approximate 
value  of  Q  required  by  the  magnitude  correction  of 
Gutenberg  (194S).  If  we  compare  the  diminution  predicted 
by  spherical  spreading  and  normal  dispersion  with  the 
diminution  predicted  by  -l.H  logA  we  sec  that  the  loss 
of  amplitude  required  in  addition  to  spreading  and 
dispersion  can  be  accounted  for  with  Q  of  about  300  at 
T  *  20  seconds*  depending  on  the  group  velocity  used  in 
(4) a  Thus,  Gutenberg's  data  gave  a  Q  less  than  that  of 
other  investigators  mentioned  above. 

Figure  4  shows  the  amplitudes  to  be  expected  for 
a  plane  wave  propagating  In  a  flat  medium  at  distances 
from  0  to  10,000  km  for  Q  values  of  100,  200»  “*00  **»<* 
1000.  From  the  evidence  above,  we  believe  that  100  and 
1000  are  the  extreme  possible  limits  to  Q  at  20  second 
periods  over  all  possible  earth  paths.  It  appears  that 
knowledge  of  Q  is  extremely  important  if  observed 
amplitudes  are  to  be  corrected  for  propagation  path. 

For  stations  10,000  km  from  the  source,  amplitudes 
could  vary  by  a  factor  of  approximately  IW0  for  the 
extreme  limits  nf  Q.  *tgaln,  this  result  must  be  altered 
by  th*  realisation  that  two  paths  from  a  source  in  the 
earth  could  not  be  so  extremely  different  in  Q  pro¬ 
perties  oves  thier  entire  length;  and  we  should 
probably  accept  factors  of  i«  or  less  as  being  more 
realistic  for  the  largest  variation  in  amplitude  due 


to  attenuation  over  two  10,000  km  paths  in  the  earth, 
l-'or  instuncc,  3.5  would  be  the  ratio  of  amplitudes  at 
10,000  km  over  a  structure  with  Q  ■  400  and  another 
with  Q  »  200  (Figure  4).  Note  that  for  distances  of 
1000  km  or  icss,  amplitudes  should  vary  by  less  than  a 
factor  of  two  even  for  the  extreme  Q  differences. 

Dispersion 

The  shape  of  the  group-veloci ty  curve  affects  the 
character  of  the  recorded  Rayleigh  wave  and  therefore 
measured  amplitudes.  This  is  most  dramatic  in  comparing 
seismograms  over  oceanic  and  continental  paths,  tfe  can 
study  the  effect  carefully  If  we  construct  synthetic 
seismograms  over  various  earth  structures  and  observe 
the  amplitudes,  he  have  already  set  down  the  various 
factors  which  couprise  the  Fourier  transform  of  a  Rayleigh 
wave  from  an  explosion  at  or  near  the  surface  of  the 
earth  --  equations  (2a),  (3a),  and  (4).  Combining  these 
we  have  fror  (l): 

If  (mi)  *  -a^l‘(w)t0(y)Au(w)k*^*r  ^^cxp(»l  (kr*3*/4)  | 
•exp|-M:r/2Q(M)ll(M))  1 1  (u) 


Replacing  k  In*  m/C,  changing  from  cylindrical  spreading 
to  spherical  spreading  of  the  earth  (using  Rfl  as  the 
earth's  radius),  and  defining  the  instrument  response  as 

l(w)  *  r.(w)c"i*tu* 


II 


where  G  (u)  is  Che  anplitude  gain  anti  »(w)  i*  the  phase 
lag,  the  transforn  of  the  recorded  signal  Is 

W(m)  «  SL  "  —  i  —  —  - 

CCwi)l,2«y2  slnl^2(r/h0)  exp|wr/2Q(u>)U(w)| 

since  up  Is  positive  on  a  seisnogran  and  we  have  adopted 
the  layering  response  of  a  Median  with  the  :  axis  posi¬ 
tive  downward,  we  shift  the  phase  an  additional  *  or, 
equivalently,  we  change  the  sign  of  the  whole  expression; 
and  to  start  the  synthetic  selsnogran  at  sene  tine  later 
than  the  origin  tine,  we  subtract  t  seconds  fron  r/C, 
where  t  Is  near  the  expected  tine  of  arrival  of  the 
wave  group  with  the  highest  velocity: 

A(w)ft|,(»)AttCw)wl/2G(M)expC-i<w|r/C(f«)-thf(M)-Sff/4}) 

*  11  1  im  —  ■ 

sInl^2Cr/R0)  e»|t|ur/2l|(w)U(u;)l 


the  recorded  Rayleigh  wave  Is  then  the  Inverse  transform 
of  (S): 

#  a°'l*(<«*)Cj|('«l^u(M)M^fcl'(Ml®*P( I Im| t*t-r/G(w) | •!(-«) -tt/t | * 

||(  t  )  *  ]  — ||w 

^  t:(M),,2Rj^2  slnl,2(r/R0)  exp|wr/2Q(M)U|w) | 

(i) 
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A  program  was  written  (I.PSW)  which  takes  as  input 
ail  the  frequency-dependent  parameters  needed  in  (I)  at 
specified  frequencies,  interpolates  them  with  straight 
lines,  and  uses  a  fast  Courier  transform  algorithm  (COOL) 
t©  construct  the  signals  represented  by  (*).  P(w)  was 
taken  to  lie  the  spectrum  of  a  step  function,  l/w,  Prom 
krun©  {lifeQ)  the  phase  velocities  for  Shield,  Mid- 
Continent,  basin  and  Mange,  Alpine,  Croat  Island  Arc, 
and  Sleep  «icean  were  used  In  synthesiclng  signals  as 
given  in  fable  I.  The  velocities  were  extrapolated  to 
five  second  periods  in  all  but  the  shield  case,  Since 
we  are  now  just  Interested  In  the  possible  effects  of 
dispersion  on  recorded  amplitudes  or  magnitude  determi¬ 
nations,  we  take  the  layer  response  (t  k1^*)  and 
attenuation  ft$)  to  be  identical  for  all  these  structures 
{knowledge  of  these  would  allow  the  amplitude  measure¬ 
ment  to  be  corrected  for  them  in  practice).  The  standard 
IMM  long-period  response  {Tab  fie  lV|  was  used  in  the  synthesl 
Values  of  ty  M  were  taken  from  Isai  and  All  IlibQ)  and 
are  given  in  Table  II  here  along  with  the  layer  response 
used,  which  happens  to  be  for  a  fiutenberg  continental 
structure.  Tor  each  new  distance,  e  was  set  to  begin 
the  signals  prior  to  the  arrival  time  of  the  groups  of 
maximum  velocity,  figure  la,  then,  shows  the  seismograms 
for  each  structure  {except  the  oceanic)  at  epicentral 
distances  from  I Oft  km  to  IQ, pun  km.  They  are  distinct 
for  the  live  structures  used.  Unexpected  high-frequency 
waves  superimposed  on  the  low-frequency  first-arriving 
groups  are  an  artifact  of  the  synthesis  technique  and 
should  be  disregarded,  figure  lb  shows  the  oceanic 


signals,  Nhich  ire  quite  different  fron  the  continental 
i4* ampler  it  i  given  distance.  figure  Sc  shows  seisno- 
gratis  f@r  the  mid-continent  structure  is  in  Figure  Sa, 
except  the  Q  values  of  Table  II  at  16  and  IS  seconds 
were  doubled;  this  was  an  attempt  to  simulate  more 
closely  J8TS  seismograms,  dominant  periods  of  which  are 
characteristically  IS  seconds  or  less  even  out  to 
teleseismic  distances. 

It  is  evident  that  the  most  precise  way  of  remov¬ 
ing  effects  of  dispersion  is  to  Fourier  transform  the 
signal  to  get  Ilf  so  that  dispersion  effects  are  con¬ 
tained  solely  in  the  phase  term  and  do  not  alter  the 
amplitudes.  If  one  does  not  have  the  time  or  machine 
capability  for  this  step,  tbe  stationary  phase  approxi¬ 
mation  to  (hi  will  indicate  the  effect  of  dispersion 
over  various  paths  on  amplitude  visually  measured  at  a 
certain  period.  Ihis  approximation  (tilth,  l$ti)  can  be 
given  as  a  function  of  distance  r  and  period  Tj 


vfl  ,r| 


l/«2| 


rj 


(Ihe  appendix  shows  the  derivation  of  this.l  'the  factor 

is  in  common  with  the  expression 
for  amplitude  of  a  spectral  component  propagating  over 
the  earth;  the  new  factor  r"1 ^  applies  to  all  periods, 
b»it  the  other  factors,  except  for  a1*,  are  dependent  on 
pe? led. 


•II 


Relative  •aptitudes  at  a  certain  period  for  the 
six  different  structrues  In  Figures  Sa  and  Sb  should  be 
given  approximately,  for  any  distance,  by: 

<1U  *> 

«  'W0U§) 

T0 

where  we  have  neglected  the  very  small  variation  (slot) 
of  6,,2(T)  In  (?)  (remember  that  the  source  spectra, 
layer  response,  and  attenuation  were  made  Identical  for 
these  six  structrues  since  we  are  not  investigating 
their  effects  here).  To  check  this,  group  velocity  dls- 
person  curves  were  first  determined  from  the  synthetic 
seismograms  at  10,000  km.  These  curves,  shown  in 
Figure  fe,  are  somewhat  sinuous  due  to  the  straight-line 
interpolation  of  the  Input  phase  velocities  by  the  syn¬ 
thesis  program.  However,  II  and  dll/dl  can  be  fairly  well 
determined  at  a  period  of  10  seconds,  and  these  values 
for  the  six  structures  are  given  In  Table  III.  Also  in 
table  III  are  the  peak-to-peak  amplitudes  measured  for 
the  six  structures  at  10,000  km  for  a  period  of  10 
seconds,  the  exact  position  of  the  10  second  period 
being  determined  by  the  arrival  time  versus  peak  number 
data,  these  amplitudes  show  a  factor  of  four  between 
extremes,  dividing  tbe  measured  amplitudes  t*y  the 
structure -dependent  factors  given  by  (S|  shew Id  result 
In  equalised  amplitudes?  however,  the  actual  results, 
given  in  Table  III  also,  show  a  factor  of  approximately 
l.y  between  the  ext  femes  of  tbe  equalised  amplitudes. 


Thus  some  equal  izat ion  has  been  accomplished,  but  this 
could  only  he  done  after  group  velocity  curves  had  been 
established  and  the  arrival  time  of  the  wave  group  of 
20  second  period  had  been  determined.  If  we  disregard 
the  oceanic  structure,  the  equalization  achieved  using 
continental  group-velocity  data  is  negligible,  since 
the  measured  raw  amplitudes  on  the  continental  structures 
varied  by  only  a  factor  of  2.2  at  the  most.  The  experi¬ 
ment  was  repeated  at  a  period  of  30  seconds  for  the 
continental  structures,  and  the  results  (not  listed) 
showed  no  benefit  In  applying  (8)  as  a  means  of 
equalizing  measured  amplitudes,  he  conclude  then  that 
amplitude  equalization  for  dispersion  effects  is 
beneficial  only  if  oceanic  paths  are  involved. 

If  we  take  from  17)  the  distance-dependent  terms 
of  sin'^fr/M^)  for  spherical  spreadlnf  and  r*1^2  for 
the  stationary  phase  approximation,  the  amplitude- 
distance  relation  for  a  given  period  Tt|  using  these  is 
shown  In  figure  7.  In  the  case  when  amplitude  is  measured 
out  an  Airy  phase  (group  velocity  minimum  or  maxinun) 
the  stationary-phase  approximation  is  invalid  and  r’^fc 
l1*  tep laced  hy  fttath,  I9t»8);  the  effect  of  this 

plus  spreading  on  the  Airy  phase  amplitude  is  also 
*h«wn  in  figure  7,  and  it  is  clear  that  at  teleseismic 
distances  use  of  the  Airy  phase  will  give  0.2  -  0.3 
higher  values  of  M  than  amplitudes  measured  at  periods 
where  the  group  velocity  curve  is  not  flat.  When  the 
absorption  term  1 1  I®)  is  added  to  the  spreading  and 
stationary-phase  effects,  the  result  I-  further  dimiuu- 
t  ion  as  shown  in  figure  ®  also,  A  ®|  of  37S»  at  I  *  In 


seconds  was  used  in  computing  the  absorption  effect; 
this  is  approximately  the  same  value  as  used  in  syn¬ 
thesizing  the  seismograms  of  Figure  5.  Also  U  was 
taken  equal  to  tt  although  in  fact  it  varies  from  2.73 
to  3.08  km/sec  for  the  five  continental  structures.  The 
line  in  Figure  7  which  represents  all  the  distance- 
dependent  factors  of  (7)  should  then  be  the  proper 
distance-correction  factor  for  use  on  the  synthetic 
seismograms  of  this  report,  or  on  any  real  siesmograms. 
We  see  that  log  A  is  not  a  linear  function  of  log  A 
over  the  entire  range  of  1°  -  100°.  But  for  the  distance 
range  A>15°,  Gutenberg's  (1945)  correction  factor  of 
-1.661ogA  fits  the  theoretical  diminution  quite  well  as 
shown  in  Figure  7,  where  the  intercept  has  been  arbi¬ 
trarily  set  for  the  Gutenberg  correction.  For  A<15° 
von  Seggern  (1970)  has  suggested  -1.091ogA  as  a  better 
fit  to  data  in  the  western  United  States,  and  this 
slope  conforms  well  to  the  theoretical  diminution  at 
A<15°  with  Q  =  375  in  Figure  7.  There  is  a  0.33  magni¬ 
tude  unit  hiatus  at  A  =  15°  between  the  Gutenberg  line 
and  the  von  Seggern  line;  the  reason  for  this  is  given 
by  von  Seggern  (1970).  The  same  theoretical  curve,  only 
with  Q  =  134,  which  is  the  value  actually  determined  in 
the  Western  United  States  by  von  Seggern,  would  still 
agree  well  with  the  -1.091ogA  correction  out  to  15° 
because  the  attenuation  factor  is  small  in  either  case 

(Figure  4).  Thus,  the  stationary  phase  approximation 
-1/2 

giving  r  amplitude  decrease  may  be  valid  to  quite 
close  distance  even  though  the  signal  is  not  yet 
dispersed.  (Figure  5  for  1000  km.)  Good  agreement  in 
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slopes  between  theoretical  diminution  and  empirical 
results  for  1°<A<100°  has  been  found.  Obviously,  the 
application  of  (8)  to  find  the  relative  amplitudes  at 
different  periods  on  one  seismogram  fails  at  close 
distances  where  the  wave  groups  of  differing  periods 
have  not  sorted  themselves  out. 

Measuring  maximum  amplitudes  from  synthetic  seis¬ 
mograms  should  enable  us  to  relate  variations  in  M  to 
the  known  composition  of  these  seismograms.  The  author 
calculated  equivalent  ground  motion  of  the  maximum 
peak-to-peak  amplitude  on  the  seismograms  of  Figure  5 
hy  dividing  the  amplitude  by  the  instrument  response 
(Table  IV)  for  the  period  of  the  maximum  excursion.  The 
amplitudes,  corrected  for  instrument  response  and 
periods  are  given  in  Table  V.  The  logarithms  of  the 
amplitudes  are  plotted  in  Figure  8  versus  the  logarithm 
of  A,  the  cpiccntral  distance  in  degrees.  The  tele- 
seismic  decay  of  amplitude  predicted  hy  Gutenberg  and 
the  regional  decay  by  von  Seggcrn  arc  also  plotted, 
along  with  the  theoretical  decay  from  Figure  7  when 
spreading  (sin"1//2[r/RQ])  ,  absorption  (cxpf-wr/2Q  (w) 

U  (oj)])  and  the  distance  dependence  of  the  stationary 

-1/2 

phase  approximation  (r  )  arc  taken  into  account. 
Gutenberg’s  linear  relation  fits  the  synthetic  data  as 
well  as  or  better  than  the  arcuate  relation  predicted 
from  theory  at  over  1600  km  distance.  In  the  stationary- 
phase  approximation  no  account  was  taken  of  the  chang¬ 
ing  period  of  the  maximum  excursion  which  is  evident 
in  the  measured  periods  of  Table  V,  Proper  use  of  the 
stationary-phase  approximation  requires  that  we  evaluate 
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the  frequency-dependent  factors  in  equation  (7)  as  well 
as  r  for  a  given  set  of  measurements  at  various 
distances  when  the  dominant  period  is  changing.  We  must 
consider  the  nearly  three-fold  increase  of  the  medium 
response  e0Aj^k  ^  (Table  II)  from  20  to  10  seconds,  the 
two-fold  decrease  of  the  source  amplitude  spectrum  P  (<*>) , 
which  we  assumed  to  be  1/w,  from  20  to  10  seconds,  and 
three-fold  increase  of  Tq"3/2  from  20  to  10  seconds.  We 
cannot  evaluate  U/(dU/dT)^2  accurately  at  10  seconds 


for  the  structures,  but  its  effect  must  be  to  offset  the 
result  of  the  variation  in  the  other  factors  because  the 


measured  amplitudes  on  the  synthetic  continental  seismo¬ 
grams  are  in  excellent  agreement  with  the  predicted 
arcuate  relation  over  the  entire  distance  range,  even 
though  this  relation  strictly  applies  to  a  constant 
dominant  period. 


The  variation  in  the  factors  mentioned  above  ex¬ 
plains  the  low  oceanic  amplitudes  in  Figure  8  because 
the  dominant  period  of  the  oceanic  seismograms  is  5  to 
10  seconds  longer  than  that  of  continental  ones.  The 
amplitude  diminution  for  the  oceanic  structure  nearly 
follows  a  -logA  relation  and  differs  considerably  from 
Gutenberg's  -1.661ogA  relation. 


Two  other  important  lines  are  included  in  Figure  8. 
One  is  the  -1.09  logA  relation  of  von  Seggern  (1970)  in 
its  correct  position  relative  to  Gutenberg's  curve. 

(As  in  Figure  7,  this  results  in  a  0.33  hiatus  at 
A  =  15°  in  the  distance  correction  factor.)  Evidently, 
the  amplitudes  from  synthetic  seismograms  of  this 
report  do  not  agree  with  those  from  NTS  seismograms. 
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The  soljrf  line  merely  confirms  this  conclusion,  because 
it  is  a  least-square  fit  to  the  NTS  data  out  to  A  ■  40° 
(von  Seggern,  1971)  and  is  a  -,901ogA  relation  (with 
arbitrary  intercept  here).  The  seismograms  of  Figure  Sc 
were  synthesized  with  higher  Q*s  at  10  to  IS  second 
periods  in  an  attempt  to  more  closely  approximate  the 
NTS  data  on  both  sides  of  the  IS0  division.  Results  of 
measuring  amplitude  on  these  are  shown  by  triangles  in 
Figure  8,  It  is  evident  that  Q's  may  need  to  be  increased 
further  and  that  the  Airy  phase  may  need  to  be  more 
distinct  before  agreement  is  attained  with  the  nearly 
-logA  relation  of  NTS  amplitudes.  Also,  changing  the 
source  time  function  from  the  step  used  here  to  a 
decaying  pulse  should  improve  the  agreement. 

The  problem  of  the  dominant  period  taking  large 
jumps,  as  for  the  shield  structure  between  200  and 
400  km  (Table  V  and  Figure  8) ,  which  causes  aberrations 
from  a  smooth  amplitude-distance  relation,  and  the 
vagaries  of  the  period  measurement  combine  to  make  the 
traditional  Mg  calculation  an  arbitrary  and  capricious 
datum  for  NTS  Rayleigh  waves,  especially  since  the  LRSM 
instrument  response  slopes  at  essentially  12  db/oct 
below  a  20  seconds  period.  A  magnitude  estimate  based 
on  the  spectrum  of  the  signal  should  be  more  reliable. 
Using  spectra  and  keeping  the  period  invariant  would 
require  an  amplitude-distance  relation  different  from 
any  of  the  three  presented  in  Figure  8;  its  form  would 
be  (from  equation  (5)): 

W  (u)  «  R"1/2  sin'1/2  (r/RQ)  exp[-ur/2Q(u)U(u) ] 


(9) 
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Ml  MU) It  AND  BOXCAR  AT  COMMON  STATIONS 


Recently  livernden  unu  Pi  Ison  (1971)  have  suggested 
that  the  difference  in  the  M^-vs-yield  behavior  for 
Amchitka  Island  and  NTS  can  he  attributed  solely  to 
the  fact  thut  measurement  of  the  surface  wave  amplitudes 
is  made  at  much  shorter  periods  for  NTS  events.  They 
state  that  at  a  given  yield  amplitude  spectra  arc  equal 
at  20  second  periods  for  the  two  sites,  hut  they  do  not 
show  individual  spectra  or  indicate  which  stations  were 
used  to  form  this  conclusion.  We  will  show  here  spectral 
computations  which  contradict  this  supposed  equivalence 
at  20  second  periods. 

Wc  have  chosen  to  compare  the  spectra  of  MII.KOW 
(Amchitka)  and  HOXCAK  (NTS)  at  five  tcleseismic  stations 
common  to  both  events:  NP-NT,  WII2YK,  PC2BC,  RK-ON,  and 
IIN-MH.  By  using  common  stations,  wc  should  he  minimizing 
effects  of  the  medium  response  factor,  c^V^k1^2.  By 
computing  spectra,  wc  should  be  eliminating  any  disper¬ 
sion  conditions  which  might  affect  visual  amplitudes. 

Two  of  the  above  stations,  NP-NT  and  WII2YK,  were  nearly 
equidistant  from  M1LKOW  and  HOXCAR  and  arc  particularly 
important  in  the  comparison. 

The  spectra  were  computed  from  the  nnfiltcrcd  seis¬ 
mograms,  digitized  at  1  sample  per  second,  using  a  fast 
I'ouricr  transform  subroutine.  Lengths  of  128,  256,  or 
512  points  (seconds)  were  used.  Since  the  tune  windows 
selected  for  spectral  computation  were  never  exactly 
one  of  these  even  powers  of  two,  zeros  were  added  to 
the  signals  and  the  computed  .amplitude  spectra  were 
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^ultipled  by  the  factor 


where 

X  *  itunber  of  point 5  in  chosen  tine  window  of 
si  glut  I 

X„  9  nunber  of  zeros  added  to  make  even  power 
of  two 

to  make  then  valid  estimates  within  the  actual  signal 
window  chosen.  The  calibrations  for  the  seismograms 
were  carefully  checked,  and  the  spectral  outputs  were 
scaled  in  millimicrons.  The  instrument  response  has 
not  been  removed,  and  the  spectra  represent  true  ground 
motion  ut  a  period  of  25  seconds  only.  One  seismogram 
has  an  uncertain  calibration  ••  MII.ROW  at  WII2VK.  This 
was  due  to  the  fact  that  the  sine-wave  calibrator  was 
inoperative  for  some  time  before  the  MII.KOW  event  and 
until  the  station  closed.  Magnification  for  this 
scimsogram  had  ti  be  figured  from  a  calibration 
obtained  about  two  weeks  prior  to  the  event. 

The  spectra  of  the  two  events  at  the  five  common 
stations  are  shown  in  figure  9,  along  with  plots  of 
the  actual  signals.  Signal  length  in  each  case  was 
controlled  by  the  group  velocities  shown  at  the  start 
aud  end  of  the  signal  plots. 
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The  spectra  Must  be  equalised  for  yield  and  dis* 
tance  differences,  and  we  will  do  this  at  the  20  second 
periods  only.  The  yields  for  BOXCAR  and  M1LR0N  are  1.2 
and  1.0  Mt,  respectively.  Since  amplitude  at  long 
poriods  is  proportional  to  yield  (von  Seggtrn  and 
l.uabert,  1969),  we  Multiply  the  BOXCAR  amplitudes  by 
the  ratio  of  1.0  over  1.2.  To  equalise  the  BOXCAR 
anplitudes  to  the  MILROtf  distance,  we  use  equation  (9) 
with  Q  •  300  and  0  -  t.  Recall  that  we  deterained  that 
this  value  of  Q  satisfied  Gutenberg's  data.  The  group 
velocity  is  arbitrarily  set  equal  to  %  in/ sec.  liquation 
(9)  applied  to  the  BOXCAR  anplitudes  neans  that  we 
Multiply  by  the  factor 


sin  (rB/Rc) 

»1»  (r,/R„) 


c*Pl-(r|,-rM)/6000| 


where  subscripts  B  and  M  refer  to  BOXCAR  and  NII.ROW, 
respectively.  The  results  of  applying  these  yield  and 
distance  corrections  arc  shown  in  Figure  9  also.  At 
NI'-NT  and  WI2YK  the  change  is  snail,  and  the  difference 
in  anplitudes  at  20  seconds  between  BOXCAR  and  MII.ROW 
is  evident.  There  i i  approximately  a  factor  of  two  at 
NI*-NT  and  a  factor  of  six  at  MI2YK.  Because  of  the 
calibration  difficulty  stated  above,  the  factor  for 
MI2YK  nay  he  erroneous.  Since  the  factor  of  six  should 
be  attributed  to  differences  in  Q  over  the  two  paths 
to  MI2YK,  a  check  with  Figure  4  at  3000  km  distance 
shows  that  only  unrealistically  disparate  Q  values 
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could  satisfy  this  anplltudc  data.  The  factor  of  two  at 
NP-NT,  however,  could  he  explained  hy  reasonable  Q 
values  of  200  for  the  MICRON  path  and  700  for  the 
BOXCAR  path.  Comparison  of  20-sccond  amplitudes  at  the 
other  stations,  PG2BC,  RK-ON,  and  IIN-MI;,  is  not  as 
simple  because  of  the  distance  differences  to  the  two 
sites.  Again,  though,  the  BOXCAR  paths  to  these  stations 
must  have  significantly  higher  Q's  to  satisfy  the  data. 

Using  common  tclcscismic  stations,  we  find  that 
the  Mg  (based  on  the  spectral  amplitude  at  20  seconds 
period)  of  MICRON  is  distinctly  lower  than  that  of 
BOXCAR  and  therefore  that  the  MR  vs  yield  relation  for 
NTS  events  cannot  be  duplicated  by  Amchitka  events, 
due  to  these  attenuation  differences. 
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CONCLUSIONS 


Hopefully,  we  have  illuminated  the  more  important 
quantities  which  cause  variation  in  Rayleigh  wave  ampli¬ 
tudes  from  site  to  site.  These  quantitites  arc  the 
medium  amplitude  response  AR,  the  attenuation  due  to 
effective  Q,  and  the  stationary-phase  approximation  for 
dispersed  waves.  All  three  deserve  consideration  when 
visual  amplitudes  are  used  to  calculate  Mg ,  but  only 
the  first  two  when  spectra  arc  used.  Another  factor  which 
may  be  important  is  multipathing ,  but  it  was  not  con¬ 
sidered  here  because  it  cannot  be  studied  in  a  general 
sense.  Capon  (1970)  shows  the  multiple  paths  inferred 
from  analysis  of  LASA  recordings  of  distant  events  along 
with  real  and  synthetic  seismograms  exhibiting  multi¬ 
ple-path  effects,  and  this  facet  of  propagation  may  be 
significant  in  relation  to  magnitude  estimation. 

It  was  found  that  the  medium  response  factor  AR 
could  vary  by  at  most  approximately  a  factor  of  three 
(0.5  on  M  scale)  between  continental  structries  and 
oceanic  structrues.  lixactly  how  well  the  amplitude  at 
a  given  site  can  be  related  to  this  factor  depends  on 
how  rapidly  the  structure  is  changing  in  the  vicinity 
of  the  site.  It  is  possible  that  oceanic  island  sta¬ 
tions  would  record  amplitudes  almost  identical  to  ocean- 
bottom  stations  nearby  on  ideal  oceanic  structure 
because  the  horizontal  dimension  of  the  transition  to 
the  island  structure  is  less  than  a  wave-length  or  two 
of  the  periods  of  interest.  We  believe  that  as  much  as 
0.3  of  the  large  negative  magnitude  differences  from 
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the  average  event  magnitude  of  island  stations  AD-iS 
(SDL  shot  report,  BRONZE)  and  IIW-IS  (SDL  shot  report, 
BILBY)  can  he  explained  by  the  medium  response  factor, 
since  these  were  the  only  oceanic  sites  for  the  two 
events . 

The  effect  of  Q  differences  can  be  very  signifi¬ 
cant,  increasing  with  telcscismic  distance.  Lor  this 
reason,  when  Q  is  unknown  for  periods  of  10  to  50 
seconds,  surface  wave  data  from  regional-distance 
stations  would  be  weighted  more  in  discrimination  work 
and  for  estimating  yields  and  seismic  moments.  This 
conclusion  was  substantiated  by  comparing  MILROW  and 
BOXCAR  spectra  at  common  stations. 

Dispersion  of  the  Rayleigh  wave  has  a  pronounced 
effect  on  amplitude  between  oceanic  and  continental 
structures  --  a  factor  of  three  or  more.  It  was  pre¬ 
viously  stated  that  AD- IS  and  IIW-IS  had  anomalously 
low  magnitudes  for  NTS  events;  they  were  actually 
about  0.8  magnitude  units  lower  than  the  event  aver¬ 
ages.  We  attributed  about  0.3  of  this  to  the  medium 
response  factor  above;  and  here  we  attribute  the 
remaining  0.5  to  dispersion  effects  such  as  illustrated 
in  figure  8  for  the  synthetic  oceanic  seismograms 
at  the  distances  in  question,  about  40°.  The  lower 
(approximately  0.4)  surface  wave  magnitudes  of  LONG 
SHOT  across  western  North  America  relative  to  more 
easterly  parts  (Lambert  ct  al.,  1969)  might  be  in  part 
attributed  to  the  fact  that  paths  to  the  Western  sta¬ 
tions  comprised  mostly  oceanic  structure.  If  only 
continental  paths  arc  under  consideration,  the 
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dispersion  effect  is  not  significant,  as  shown  by  the 
clustering  of  magnitude  values  in  Figure  8  for  tele- 
seismic  distances  at  least.  Dispersion  corrections 
bused  on  the  stationary-phase  approximation  for  con¬ 
tinental  paths  probably  docs  not  warrant  the  effort, 
especially  since  this  approximation  is  difficult  to 
evaluate  in  practice. 

The  use  of  spectra  to  estimate  M  would  obviate  the 

•» 

requirement  for  a  dispersion  correction  for  every  path 
and  would  also  allow  one  to  escape  pitfalls  such  as 
sudden  jumps  in  the  predominant  period  with  distance. 
These  jumps  can  cause  magnitudes  at  regional  distances 
to  vary  erratically  if  one  adheres  to  the  rule  of 
measuring  the  maximum  amplitude  on  the  seismogram.  This 
pitfall  is  somewhat  analogous  to  the  problem  with  body- 
wave  magnitude  when  different  branches  of  the  travel 
time  curve  lie  close  together.  Also  for  this  reason, 
visual  surface  wave  data  at  regional  distances  should 
be  weighted  less  than  that  at  tclcseismic  distances; 
this  balances  the  ill  effects  of  Q  on  tclcseismic 
surface-wave  amplitudes.  If  spectra  arc  used,  however, 
the  regional  data  should  still  be  more  reliable. 

Finally,  we  wish  to  emphasize  the  advantages  of 
using  longer  periods,  say  SO  seconds,  in  determining 
Ms .  We  already  know  that  due  to  differences  in  source 
parameters  (depth,  dimensions,  time  function)  the  M,.  vs 
m^  criterion  seems  to  discriminate  even  better  at  longer 
periods  than  at  20  seconds  (Molnar  ct  al.,  1009).  In 
addition,  all  the  propagation  factors  which  affect  M 
values  as  discussed  in  this  report  arc  much  less 
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vsttNMt  anottfi  typical  earth  iirwiur«i  at  periods  tf 
SO  seconds  than  at  periods  of  2ti  seconds.  The  amplitude 
response  of  a  layered  medium,  Ag,  differs  hy  only  a 
factor  of  at  most  I.S  at  SO  seconds  periods  for  the 
same  group  of  structures  for  which  Figure  I  showed  a 
facto*  of  almost  S.tt  at  20  siconds  periods.  Available 
evidence  suggests  that  the  effective  Q  for  SO  seconU 
waves  lies  between  100  and  200  globally,  while  this 
quantity  for  20  second  waves  may  be  anywhere  from  100 
to  over  1000  because  of  the  extremely  heterogeneous 
nature  of  the  crust.  The  character  of  dispersion  at 
SO  seconds  is  less  variable  than  at  20  seconds,  but 
equalization  procedures  should  still  be  applied  If 
oceanic  and  continental  paths  are  involved.  Again  this 
could  he  done  by  calculating  the  stationary-phase- 
approxination  factor  or  by  computing  spectra  of  the 
signals.  Thus,  at  longer  periods  we  can  take  advantage 
of  more  consistent  behavior  of  surface  waves  over  the 
large  range  of  possible  earth  structures,  thereby 
reducing  uncertainties  in  and  yield.  Also,  use  of 
longer  periods  which  are  the  first-arriving  wave  groups 
nay  reduce  concern  over  multiple-path  effects  on  ampli¬ 
tudes.  It  has  yet  to  be  shown  that  the  signal-to-noise 
ratio  for  Sfl-second  waves  from  explosions  can  be  made 
to  exceed  that  for  2l»-seeond  waves  though. 
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Phase  Velocities  Taken  From  Brune  (1969)  and  Used  in 
the  Seismogram  Synthesis  of  this  Report 
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TABLL  II 


Values  of  Q  (Tsai  and  Aki,  1969)  and 
Response  for  a  Gutenberg  Continental 
Menahem  and  llarkrider,  1964)  used  in 

thesis  of  this 


Values  of  the  Medium 
Barth  Model  (Ben- 
the  Seismogram  Syn- 
Report 
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j-qua lined  amplitude  is  equal  to  the  measured  amplitude  divided  by  U(dU/dT) 


TABLE  IV 


Relative  System  Magnification  for  LRSM  Long-Period  Recording 


Relative 

Period  (sec)  System  Magnification 


8 

.100 

9 

.138 

10 

.188 

11 

.245 

12 

.310 

13 

.390 

14 

.480 

15 

.580 

16 

.690 

17 

.790 

18 

.880 

19 

.920 

20 

.960 

21 

.970 

22 

.980 

23 

.990 

24 

.995 

25 

1.000 

26 

.990 

27 

.960 

28 

.920 
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!•  igurc  1.  Relation  ol  the  Ray  1  c  i  gh  -  wave  amplitude*  response 
factor  to  the  product  of  group  and  phase  velocity  ( twenty - 
second  period  of  the  fundamental  mode)  for  various  earth’ 
structures . 
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Figure  2.  Relation  of  the  Rayleigh-wave  amplitude  response 
factor  to  group  velocity  (twenty-second  period  of  the 
fundamental  mode)  for  the  same  earth  structures  as  in 
Figure  1. 
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I'igure  3,  Relation  of  the  Rayleigh-wave  amplitude  response 
factor  to  phase  velocity  (twenty-second  period  of  the 
fundamental  mode)  for  the  same  earth  structures  as  in 
Figure  1. 


I'igurc  4.  Amplitude-distance  relation  of  surface  waves  for 
various  values  of  Q. 
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figure  5a.  Synthetic  fundamental-mode  Rayleigh-wave 
signals  for  five  different  continental  structures  at 
various  cpicentral  distances. 
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signals  for  five  different  continental  structures  at 
epicentral  distances. 


Figure  5a  (cont'd.).  Synthetic  f 
signals  for  five  different  conti 
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Figure  5a  (cont'd.).  Synthetic  fundamental -mode  Rayleigh- 
signals  for  five  different  continental  structures  at  vari 
epicentral  distances. 


different  continental  structures 


MAXIMUM 

AMPLITUDE 

43.3 


4  7 


Figure  5a  (cont'd.).  S 
signals  for  five  diffe 
eDicentral  distances. 
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Figure  5a  (cont'd.).  Synthetic  fundamental -mode  Rayleigh-wav 
signals  for  five  different  continental  structures  at  various 
epicentral  distances. 
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Figure  5a  (cont’d.).  Synthetic  fundamental -mode  Rayleigh- 
signals  for  five  different  continental  structures  at  varx 
epicentral  distances. 
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Figure  5a  (cont'd.).  Synthetic  fundamental -mode  Rayleigh-wave 
signals  for  five  different  continental  structures  at  various 
epicentral  distances. 


MAXIMUM 

AMPLITUDE 


3 


Figure  Sa  (cont'd.).  Synthetic  fundamental -mode  Rayleigh-wav 
signals  for  five  different  continental  structures  at  various 
epicentral  distances. 
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Figure  5a  (cont*d.).  Synthetic  fundamental -mode  Rayleigh- 
signals  for  five  different  continental  structures  at  vari 
epicentral  distances. 
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Figure  5b  (cont'd.)-  Synthetic  fundamental -mode  Rayleigh  wave 
signals  for  a  deep  ocean  structure  at  various  epicentral 
distances. 
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Figure  5c.  Synthetic  fundamental-mode  Rayleigh-wave 
signals  for  a  high-Q  mid-continental  structure  at  various 
cpiccntral  distances. 
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Figure  5c  (cont'd.).  Synthetic  fundamental -mode  Rayleigh-wave 
signals  for  a  high-Q  mid-continental  structure  at  various 
epiccntral  distances. 


MAXIMUM 


Figure  5c  (cont|d.)  Synthetic  fundamental-mode  Rayleigh-wave 
signals  for  a  high-Q  mid-continental  structure  at  various 
epicentral  distances. 
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Figure  6.  Group  velocity  versus  period  as  measured  on  the 
10,000  km  synthetic  seismograms  of  Figures  5a  and  5b. 
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ure  6  (cont'd.).  Group  velocity  versus  period  as  measured 
the  10,000  km  synthetic  seismograms  of  Figures  5a  and  5b. 
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Figure  6  (cont'd.)*  Group  velocity  versus  period  as  measured 
on  the  10,000  km  synthetic  seismograms  of  Figures  Sa  and  Sb. 
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Figure  6  (cont'd.).  Group  velocity  versus  period  as  measured 
on  the  10,000  km  synthetic  seismograms  of  Figures  5a  and  5b. 


Figure  6  (cont’d.).  Group  velocity  versus  period  as  measured 
on  the  10,000  km  synthetic  seismograms  of  Figures  5a  and  5b. 


Figure  8.  Measured  maximum  amplitudes  from  the  synthetic 
seismograms  of  Figures  5a,  5b  and  5c  compared  to  theoretical 
and  empirical  amplitude-distance  relations. 
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I'igurc  lJ.  Spectra  of  BOXCAR  and  MILKOW  Rayleigh  waves  at 
NI’-NT,  WII2YK,  PG2UC ,  KK-ON,  and  IIN-Mli. 
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Figure  9  (cont'd.).  Spectra  of  BOXCAR  and  MIL ROW  Rayleich 
at  NP- NT ,  WII2YK ,  FC2BC ,  RK-ON,  and  IIN-MI-. 
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1-iKure  !>  (contM.).  Spccti;i  of  BOXCAR  ;md  MII.ROK  Ravloich  iv.ivos 
:‘t  Nl*  NT,  IVII2VK ,  I'ddBC,  RK-0\,  and  IIN-MH.  M 
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DISCRETE  AMPLITUDE  SPECTRUM,  mfj 


APPENDIX 


DERIVATION  OF  STATIONARY-PHASE  APPROXIMATION 


Equation  (6)  in  the  text  gave  the  expression  for 
the  fundamental-mode  Rayleigh-wave  signal  from  an  explo¬ 
sion.  Let  du>  =  Udk  and  w/C  *  k  in  that  expression  to 
obtain: 


w(r,t) 


sin1/z(r/R0) 


PeoARGUu1/2exp[i(a)t-kr)-<J.-5ir/4] 
C1/Z  exp[-u>r/2Q(u>)U(ui)] 


This  can  be  compared  with  a  similar  expression  in  Bath 
(1969,  p.  44,  equation  (9)).  We  must  now  assume  that  all 
the  factors  in  the  integrand  except  the  eikonal  function 
vary  slowly  with  wavenumber  k.  The  stationary-phase 
approximation  for  (Al)  is  then  given  by  reference  to 

o 

Bath  again  (p.  45,  equation  (13)): 


w(k,r) 


a  3Pe0ARGUa) 2cos  (oot  -  k  r- 4.  -  5tt/ 4+ir/ 4 ) 


C1/2R*/2  sin1/2(r/RQ)  exp(-wr/2QU) 


Ttr  d2oj 
2U  dk2 


1/2 


where  we  have  replaced  t  by  r/U  in  Bath's  expression. 

The  positive  sign  applies  to  normal  branches  of  the 
dispersion  curve,  the  negative  sign  to  inverse  branches. 
We  desire  to  express  u>  in  terms  of  period  TQ  rather 
than  wavenumber  kQ.  Since  U  =  dio/dk,  we  can  write: 


or 


. 2 

d  u) 

- 7 

dx 


U 


dU 


From  the  definition  w  =  2it/T 


du> 


27rdT 

—’T 

T 


Thus  (A3)  becomes 


(A3) 


UT2  dU 

TfT"  HT 


(A4) 


Substituting  (A4)  in  (A2)  and  using  k  =  2tt/CT  and  wt  *  27rr/UT: 


a3Pe  AdGUu)1/2  cos(27rr/UT-2TTr/CT-<t>-S7r/4+Tr/4) 
o  K  "  -  - 

C1/2Rq1/2  sin1/2(r/RQ)  exp(-Ttr/ QUT) 


Simplifying,  we  have 

7r1/2a3PeoARGU  cos  ( (27rr/T)  (1/U-l/C) -$-<;it/4+it/4] 

1  21/2C1/2Rq1/2  sin1/2(r/RQ)  exp(-irr/QUT)  T3/2  |jjl 

The  form  of  equation  (8)  in  the  text  follows  if  we  assume 
that  all  the  factors  are  constant  except  those  determining 
the  dispersion,  namely,  U  and  C. 
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